Coherent atomic soliton molecules for matter-wave switching 
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We discuss the dynamics of interacting dark-bright two-dimensional vector solitons in multicom- 
ponent immiscible bulk Bose-Einstein condensates. We describe matter-wave molecules without a 
scalar counterpart that can be seen as bound states of vector objects. We also analyze the possibility 
of using these structures as building blocks for the design of matter-wave switchers. 
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Introduction. Solitons are robust wave-packets able to 
maintain their shape when propagating in different me- 
dia and under mutual collisions. The existence of such 
clastically interacting localized waves is an essential prop- 
erty of integrable nonlinear equations that have an infi- 
nite number of conservation laws Integrable two- 
dimensional equations, such as Kadomtsev-Petviashvili, 
Davcy-Stcwardson, Zakharov-Manakov, sine-Gordon and 
others have two-dimensional solitons with complex inter- 
actions. Some nonintegrablc systems may have solitary 
waves - localized coherent structures with almost clastic 
interactions. An interesting open question is the con- 
struction of complex molecule-like coherent matter-wave 
structures, i.e. superpositions of solitons leading to sta- 
ble bound states with molecule-like behavior. 

In this paper we present nontrivial nonlinear phe- 
nomena in multicomponent Bose-Einstein condensates 
(BECs) described by coupled two-dimensional nonlinear 
Schrodingcr equations (NLSEs). The balance of disper- 
sion and nonlinear interactions in BECs leads to differ- 
ent types of nonlinear coherent excitations (see e.g. the 
experimental papers [HI] or the reviews 0, El). Soli- 
ton molecules have been considered in the propagation of 
optical beams in nonlinear media with saturable nonlin- 
caritics 

(UHll, it being very difficult to construct even 
metastablc long-living soliton clusters with local interac- 
tions such as those present in ordinary BECs. 

Here we will construct soliton molecules using soliton- 
bubble bound states (i.e. two-dimensional extensions of 
the dark-bright soliton pair) as bricks to construct matter 
wave aggregates. These can also be viewed as a combina- 
tion of the extension of a Kadomtsev-Petviashvili soliton 
[TB1 | and a bright NLS soliton or as a vortex pair super- 
imposed with two density peaks. 

Matter-wave trains with a finite number of one- 
dimensional bright solitons are stable due to the pres- 
ence of the trap [HI. However, the idea does not work 
for higher dimensions due to the blow-up phenomenon 
Q. With defocusing nonlinearitics, dark solitons always 
repel each other and cannot form bound states fl7j |. 

We will show how multicomponent homonuclear BECs 



in the immiscible regime allow for the construction of 
robust novel types of solitonic molecules. These matter- 
wave clusters display phase-dependent properties due to 
their coherent nature and can be used for constructing 
nonlinear matter- wave switchers. 

Physical system and model equations. We will consider 
two-component BECs with atoms in two hyperfme states 
1) and 1 2) in the immiscible regime and consider droplets 
of atoms in component |2) to be phase separated from a 
component |1) assumed to have a much larger number 
of particles. When tightly confined along one direction 
these systems are ruled in the mean field limit by 



' dt 



9jk\4>k\ 2 -fi + Sj\ ^j, (i) 



ifc=l,2 



for j = 1,2. Without loss of generality we work in di- 
mensionless units, with chemical potential fj, = 1, and 
Si = 0. Immiscibility implies that gf 2 > 9\\9ii- The 
normalization for ip 2 is given by J R \ip2\ 2 = 2(022/00)^2, 
where 022 and N 2 are the s-wave scattering length and 
number of atoms in |2). Finally oq — ^Jh/mui^ is the 
length-scale in which spatial units are measured. 

Soliton molecules. In one-dimensional one-component 
NLS systems, the repulsive nature of the interaction be- 
tween dark solitons prevents them from generating bound 
states. The interactions between bright solitons in the 
absence of external effects (such as external confinement) 
depend on the phase differences, A<p = |0i — going 
from attractive for Acj> = to repulsive for Ad) = tt. A 
critical intermediate regime for A<f> exists in which un- 
stable bound states can be constructed. 

A different possibility is constructing a vector object 
including a dark soliton in one component and a bright 
soliton in another component. This dark-bright pair 
could lead to a stable bound state with a second vec- 
tor soliton of the same type when the repulsive interac- 
tions between the dark components are balanced by the 
attractive interactions between the "droplet-like" bright 
components [3|. However, when passing to higher di- 
mensions the phenomenology changes essentially due to 
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FIG. 1: (Color online) (a) Energy (E — E gs )- impulse (p = 
pi +P2) dispersion curve of solitary wave solutions of Eqs. J2]| 
for a = 1.2 and various choices of iVV The upper solid line 
corresponds to A^2 = 0, i.e. the JR dispersion curve, shown for 
comparison. To ease comparison we subtract the ground state 
energy E gs in the plot. The dashed black line and the dashes- 
dotted line correspond to iV2 = 8 and N2 = 20, respectively, 
(b) Velocity as a function of momentum for N2 = 20; the 
red (monotonically decreasing) curve corresponds to the JR 
case; the blue (passing through the origin) curve correspond 
to nonzero N2. 



the fact that the transverse instability of the dark soliton 
leads to the formation of vortex pairs of opposite circu- 
lation. Therefore, the more natural building blocks for 
a bound state are vortices of the first component host- 
ing a "droplet" of the second component. In the scalar 
case, all movin g tw o-dimensional coherent states were 
found in Rcf. [19}. These correspond to vortex pairs 
of opposite circulation and rarefaction pulses. A similar 
phenomenology arises in the two-component case in the 
miscible regime [2(| • As the velocity of the solitary wave 
increases, the distance between vortices of opposite circu- 
lation decreases to zero. The solutions at even higher ve- 
locity are localized density perturbations without zeroes. 
In two dimensions the sequence of solutions terminates 
with solutions approaching zero energy and momentum 
as the velocity U approaches the speed of sound. 

We will construct solitary waves with velocity U 
along the x— direction in two-dimensional two-component 
BECs in the phase-separation regime as solutions of Eqs. 
([T|) in the frame moving with the disturbance: 

iU^ = IvVx + Cl-M 2 -*! 2 )^, (2a) 



iU 



dx 
dx 



1 



VV2 + (A - a|^i| 2 - IV2I 2 ) ^2 (2b) 



together with the boundary conditions — > 1, ip2 —> 
0, as |x| — > 00. In the phase separation regime a = 
312/511 = 912/922 > 1- Here A = /U2//X1 where ^1 
and \X2 are the dimensional chemical potentials of 
and ip2- We solve numerically the discretized version 
of Eqs. @ by a Newton- Raphson algorithm combined 
with a secant algorithm to find A for a given constraint 
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FIG. 2: (Color online) Density plots of the first (left) and sec- 
ond (right) components of a stationary bound state solution 
for U = 0.1, -/V2 = 80. White color indicates maximum den- 
sity. Streamlines of the first component are shown on the left 
panel. The spatial region shown is 60 x 40 healing lengths. 



on N2 — J \tp2\ 2 dx dy . We obtain a family of solutions 
characterized by the velocity of propagation U, energy, 
E and impulse p = (pi +P2, 0), given by 



E = 



^|V^i| 2 + i|VV; 2 | 2 +c#i| 2 |^2| 2 



+ ^(1 - K'll 2 ) 2 + - M^fdxdy 
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dxdy 



,(3a) 
(3b) 



with j = 1,2. The resulting families of solutions are 
plotted in Fig. [T] for various choices of N2 together with 
the Jones-Roberts (JR) dispersion relations [19( for one- 
component condensates. For a given speed U, solitary 
solutions with higher a have lower energy and higher 
impulse. 

In contrast with the JR solutions, there is a stationary 
solitary wave with nonzero energy E gs corresponding to 
the ground state of the system with all the mass of the 
second component forming a radially symmetric "bub- 
ble" in the center of the depleted first component. As the 
velocity increases from zero, the bubble becomes oblate 
in the direction of the motion with the velocity field of 
the first component being that of a dipole (see Fig. [5]). 

There is a point on the dispersion curve where the ve- 
locity reaches its maximum - the inflection point. As 
energy and momentum increase, the velocity decreases 
and the solutions become pairs of vortices of opposite 
circulation in the first component with the second com- 
ponent filling the vortex cores. In general, bubble-like 
solutions for small E can be seen as a bound state of a 
JR rarefaction pulse and a mobile "filling" of the second 
component. Fig. W[h) shows that there is a maximum 
velocity for the propagation of these solutions (different 
from the sound speed). This is a signature of the mass 
of the second component, the heavier being the second 
component the smaller being this velocity. 

We have simulated the evolution of two bubbles set 
on a colliding course. Initially they are separated by a 
large distance, so that individually they are accurately 
represented by the solutions we found. Several possi- 
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FIG. 3: Density snapshots of the second component in co- 
herent bubble-droplet pairs collisions, before (left column) 
and after (right column) the collision, (a-b) Identical inci- 
dent bubbles with U = 0.1, N 2 = 40, p = 6.36, E = 2.99 form 
a bound state, (c-d) Incident bubbles with the same speed 
Ui = U% = 0.2 but different sizes emulate an elastic colli- 
sion (left bubble: N 2 = 40, E = 5.6, p = 14.9, right bubble: 
N 2 = 50, E = 7.7, p = 20.1). (e-f) Identical bubbles (param- 
eters being U = 0.2, N 2 = 40, E = 5.6, p = 14.9) emulate the 
collision of two pairs of well-separated vortices, (g-h) Identi- 
cal bubbles (same as those in third row) collide with an offset 
of 8 dimensionless units and form a bound state. The spatial 
region spanned is £ £ [—50, 50], y 6 [—30, 30]. 



ble outcomes of such collisions are summarized in Fig. 
[31 Almost identical slow colliding bubbles may form a 
bound stationary state even when they collide with an 
offset. Bubbles moving with large velocities may scatter 
at 7r/2 angle resembling the collision of two pairs of well- 
separated vortices of opposite circulation. Almost elastic 
collisions between these structures were observed when 
the velocities or masses of the bubbles were very different. 
A bound state is more likely to be formed when bubbles 
have similar phases of the second component and move 
slowly. In such collisions, a small fraction of the mass 
is emitted as sound waves. The outgoing bubbles are 
solitary waves as verified by energy-impulse calculations. 

Our simulations of Eqs. (TTJ) were done using two high 
order schemes: a fourth-order in space finite differences 
with a fourth-order Runge-Kutta scheme in time (FD) 
and a Fourier pseudospectral method (FS) [2l|, with 
Ax = Ay = 0.25 and At = 0.01 (FD), 0.005 (FS). 

Dynamic molecules. These solitonic molecules may 
have complex oscillatory internal modes: the states 
formed by vortices with a bright filling, exchange parts 
of the mass leading to a periodic beating, as we will show 
simulating the evolution of initial conditions: 
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(4b) 



Being r± = y]{x ± f ) 2 + y\ ± = atan (y/ (x±f)), 
Wd and Wb, respectively, the widths of the vortex cores 
and the "bright" component of these multidimensional 
dark-bright solitons and D the separation between the 
centers of the "dark" parts. The choice of Gaussian 
and hyperbolic tangent profiles for the bright compo- 
nent and vortices, respectively, is more realistic from the 
experimental point of view than numerically calculated 
eigenstates. We consider N 2 = 5 as number of parti- 
cles of the bright component. We also impose Wb = 1.0 
and that the diameter of the bright part fits exactly 
within the vortex size at one half of the maximum ampli- 
tude, which yields Wd = 2.3, calculated by the condition: 
IV>i (^^nr 2 ) I 2 = HV'iMI 2 . We have chosen gij = -1, 
with j = 1,2, to illustrate that the results are valid 
even in the boundary of the immiscibility region. In this 
configuration the two nonlinear structures move along 
a rectilinear path, keeping the distance D = 4.5 fixed. 
The bright component of the pair first tunnels to the 
empty vortex core and then beats periodically between 
the two cores similarly to the bright soliton oscillation 
in a double- well potential (see Fig. S] (top)). The vor- 
tex cores significantly change their size along the whole 
dynamics, the effective potential wells being dynamically 
modified along the direction of propagation owing to the 
inter/intra-component interactions (see the inset in Fig. 
3]). Initially (point (a)), the vortex is empty, and there- 
fore features a minimum width. As the bright component 
starts to fill up the vortex core, it becomes broader ow- 
ing to the repulsive interaction between different atomic 
species, reaching a maximum width for t = 37.5 (point 
(b)). When the bright component in a vortex decreases, 
it narrows again (see point (c) for t = 75). 

The bottom plot in Fig. 2] shows the variation of the 
maximum of the bright population within the initially- 
empty core N max (normalized to the total number of par- 
ticles in the bright component) as a function of D. As it 
can be appreciated in the graph, when the cores are close 
enough, most particles oscillate between both cores. A 
small fraction of particles remains in the empty vortex, 
because of the overlapping between modes in both cores. 
As D increases over a certain threshold the tunneling be- 
tween cores is suppressed due to the nonlincarity-induccd 
asymmetry of the effective potentials and the value of 
Nmax quickly decreases to zero. Thus, this switching 
curve, having an intermediate regime of partial tunnel- 
ing in which is possible to control N max efficiently, allows 
for the design of a "matter-wave switcher". It is also 
possible to control both N max and the beating period 
by slightly tuning the interespecies coupling coefficients 
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FIG. 4: Top: Propagation of a vortex pair of opposite circu- 
lation in component |1) with a bright soliton in component 
1 2) nested into one of the cores. The pseudocolor plot shows 
the density of the bright component. The snapshots are taken 
in the center-of-mass system with a separation of 6.25 time 
units, each showing a spatial window of 23.5 x 9.5 units of 
a larger simulation region (100x100). The initial separation 
between both cores is D = 4.5. Bottom: Switching curve 
N m ax vs. D. Inset: Evolution of Wd for the lower vortex core 
in the simulation displayed in the top. 




FIG. 5: (Color online) Density isosurface plot of the beating 
process of the bright field hosted by two equally-charged vor- 
tices. Dark cores rotate keeping the distance to its center of 
mass unaltered. A periodic beating of the bright component 
(blue surface) is observed in between both cores, as it can be 
inferred from the white contours within the inner slices. In 
this plot, D = 3.5, the grid size is 100 x 100, the spatial range 
is x, y € [—10, 10], a = 1 and the time interval t € [0, 150]. 

(gi2, <?2i), making the switcher completely reconfigurable. 
Control of g\2 and <?2i has been demonstrated previously 
in experiments (22j . The beating persists for times longer 
than the condensate lifetime, indicating that the nonlin- 
ear matter-wave switching process is very robust. 

A second interesting configuration corresponds to a 
pair of vortices of equal charges and the same initial con- 
ditions described previously in the discussion of the re- 
sults shown in Fig. 01 For this system, the vortices orbit 
following circular trajectories while the bright component 
tunnels between the vortex cores and thus follows a spi- 
ralling trajectory (see Fig. [5]). This beating is very stable 
and does not affect the circular vortex trajectories. 



Conclusions. We have discussed the dynamics of in- 
teracting bubble-droplet pairs in quasi-two dimensional 
immiscible BECs. We have found novel types of robust 
bubble-like solitons without a scalar counterpart that can 
be used to construct coherent atomic soliton molecules. 
The study of the dynamics of these stable objects has 
revealed the possibility of constructing nonlinear matter- 
wave switchers. Our ideas could be tested in future ex- 
periments. 
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